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ABSTRACT 

Solar flare accelerated electron beams propagating away from the Sun can interact with the turbulent 
interplanetary media, producing plasma waves and type III radio emission. These electron beams are 
detected near the Earth with a double power-law energy spectrum. We simulate electron beam 
propagation from the Sun to the Earth in the weak turbulent regime taking into account the self- 
consistent generation of plasma waves and subsequent wave interaction with density fluctuations 
from low frequency MHD turbulence. The rate at which plasma waves are induced by an unstable 
electron beam is reduced by background density fluctuations, most acutely when fluctuations have 
large amplitudes or small wavelengths. This suppression of plasma waves alters the wave distrubtion 
which changes the electron beam transport. Assuming a 5/3 Kolmogorov-typc power density spectrum 
of fluctuations often observed near the Earth, we investigate the corresponding energy spectrum of 
the electron beam after it has propagated 1 AU. We find a direct correlation between the spectrum of 
the double power-law below the break energy and the turbulent intensity of the background plasma. 
For an initial spectral index of 3.5, we find a range of spectra below the break energy between 1.6-2.1, 
with higher levels of turbulence corresponding to higher spectral indices. 
Subject headings: Sumflarcs - Sun: X-rays, gamma rays - Sumactivity -Sun: particle emission 



1. INTRODUCTION 

Solar flare impulsive electron events present an al- 
ternative to the more traditional hard X-ray diagnos- 
tics of poorly understood acceleration and transport 
of solar energetic electrons. While hard X-ray obser- 
vations provide insight into energ etic electrons in the 
lower dense solar atmosphere (e.g. lArnoldv et all 119681 : 



theor e tical (iZheleznvakov fc Zaitsevl I1970L iZaitsev et al.l 
119721: iMel 'Nik 19 95|) and numerical investigations 
(iMagelssen fc SmitIilll9T7l:lGrognardlll982l:lKontar' et all 
19981: lYoon et al.l 12 000: K ontarl 12001b!: iLi et al.l 12006a!: 



iDennis fc Schwartzlll989ilBrown fc KontarlteOOlTimpul 
sive s olar electron events (e.g. iLinl 119851 : iKrucker et all 
l2007f i provide crucial information about escaping elec- 
trons from the acceleration region. Because of the rather 
limited spatial resolution of past and current hard X- 
ray observations, even the spatially resolved hard X- 
ray spectrum of en ergetic electrons with RHESSI (e.g. 
lEmslie et al.l I2003D is a convolution of transport ef- 
fects and possibly electron acceleration dBrown et al.l 
l2009f ). Hence, de-convolution of the electron acceler- 
ated spectrum and acce lerator properties from hard X- 
rays (|Brown et al.ll200"6[) is a non-trivial task using both 



forwa rd-modelling dHolman et al . 2003; Kasp arova et al 



2005 ) and model-indepen dent techniques (jPiana et al 



20031 iKontar et al.1 120051 ). Solar electron impulsive 



events propagate outward through the almost colli- 
sionless pl asma of the solar corona and solar wind 
(|Linl [1985). Even with this collisionless regime the 
energetic electrons can interact with plasma via gen- 
eration and absorption of el ectrostatic plasma waves. 
In t he standard scenario (jGinzburg fc Zhelezniakovl 
119581 ) , the non- linear interaction of beam-driven plasma 
waves leads to the appearance of rather strong radio 
emission - type III solar/interplanetary radio bursts. 
The observations of type III so l ar bursts and en- 
ergetic particles dLin et all 119811: lErgun et all 119981: 
iGosling et all 120031 : IKrucker et all |2007D as well as 
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Ledenev et al.l 120041: iKrasnoselskikh et all 120071 ) provide 
strong support to the standard type III model. 

The plasma of the solar corona and the solar wind is 
a non-uniform turbulent medium with density pertur- 
bations at various length scales. The structure of the 
solar wind density fluctuations have been anal ysed us- 
ing scintillations o f small-size radio sources (e.g. lHollwegl 
119701 : lYoung|[r971l ) . In-situ measurements have also been 
used to determine the density spectrum n ear the Earth 
and between 0.3 and 1 AU with Helios (jMarsch h Tul 
119901) . While the detailed structure of the density tur- 
bulence in the inner heliosphcrc is not well established, 
the density fluctuation spectrum near the Earth seems 
close to a power-law spectrum with spectral index about 
5/3, sim ilar to earlier observations. It has been rec- 
ognized (|Rvutovl 119691 : iKarpman fc Istominl 119741 ) that 
beam-driven Langmuir waves can be effectively altered 
by even weak density gradients. Therefore the effect 
of density fluctuations on beam-driven plasma waves 
responsible for type III radio bursts has been consid- 
ered both numerically and analytically (jMelrose et all 
1986: Robinson et al. 1992: Kontar 2001a). Density fluc- 
tuations are believed to suppress plasma wav e growth 
(jSmith fc Simd 119791: iMuschietti et al.l 119851) and be 
responsible for the clumpy plas ma wave distribution 
observed in-sit u near the Earth (jGurnett et all 119781 : 
iLin et ail 119811 ). The fluctuations, whilst changing the 
distribution of plasma waves significantly, have a rather 
weak modula tion effect on th e instantane ous distribution 
of elec trons ([Kontarl l2001bD . Recently, IKontar fc Reidl 
(|2009f) have shown that the electron beam plasma in- 
teraction via Langmuir waves in the non-uniform solar 
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corona leads to the appearance of a break energy in the 
observed spectrum at the Earth and can explain the ob- 
served apparent early injection of low-energy electrons. 
However, the net effect of density fluctuations in the so- 
lar wind on the electron spectrum detected near 1 AU 
has not been addressed before. 

In this paper, we investigate the effects of background 
plasma density fluctuation on the generation and absorp- 
tion of plasma waves from a high energy solar electron 
beam travelling from the Sun to the Earth. We demon- 
strate the dependence of plasma waves on the level of 
density fluctuations, with high levels damping plasma 
waves too much to be in accordance with detected type 
III radio emission. We also show how the level of density 
fluctuations has a direct effect on the spectral character- 
istics of the electron beam near the Earth. 

2. ELECTRON BEAM TRANSPORT MODEL 

There is a variety of different processes which affect 
the propagation of high energy electr ons from the so- 
lar corona through the heliosphere (see lMerroselll99(]t as 
a review). This work focusses on the role of electron 
beam-driven electrostatic turbulence in the propagation 
and spectral evolution of energetic particles. The electro- 
static turbulence plays the dominant role for deca-keV 
electrons. The solar magnetic field expanding into the 
heliosphere quickly decreases with distance and provides 
adiabatic focussing for energetic electrons which ensures 
one dimensional (along expanding magnetic field lines) 
electron transport. To describe self-consistently resonant 
interaction of the electron distribution function f(v, r, t) 
(the number density of energetic electrons is 7i& = J fdv) 
and the spectral energy density of electron plasma waves 
W(v, r, t) (the energy density of plasma waves is / Wdk 
ergs cm -3 ) in the radially expanding magnetic field of 
the heliosphere, one can use the following equations of 
weak turbulence theory 
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The first terms at the right hand sides of Equations (|H2p 
describe the resonant interaction, u ! ve = kv of electrons 
and p l asma waves first deriv ed by iDrummond fe Pinesl 
(|1962MVedenov et al.l (|1962f ). The dispersion relation of 
plasma waves is wl(/c) = uj pe + "iv\ e k 2 / '(2w pe ), so the 
group velocity of plasma waves is duj^/dk = 3v 2 -, e /v 

in Equation d2| where VTe = \/kBT e /m e . Following 
IZheleznvakov fe Zaitsevl (|1970D ; ITakakura fe Shibahasiul 
(| 19761 ) we include collisional losses both for electrons 
and Langmuir waves. The last term of Equation 
([lj accounts for electron collisiona l Coulomb lo sses in 
fully ionized hydrogen plasma (e.g. lEmslid 119781) . 7 C = 
7rn e e 4 In A/(??igW 3 n e ) is the collisional damping rate of 



is the Landau damping of Langmuir waves by back- 
ground plasma. The last term in Equation ([2]) is 
the spontaneous wave g enera t i on, which is similar to 
IZheleznvakov fc Zaitsevl (119701) : ITakakura fc Shibahashl 
(|1976D : lHannah et all (120091) but different from the terms 
used in iLi et all (j2006bl ). We note that for large ve- 
locities (v > t>Te\/21nA) the energy loss of an electron 
to spontaneousl y generate Langmuir waves adopted by 
ILi et al.l (|2006bf ) is greater than the electron collisional 
Coulomb losses in fully ionized hydrogen plasma (last 
term of Equation [T]) . 

The second term on the left hand side of Equation (fT]) 
models magnetic field expansion from the solar corona 
into interplanetary space and the 'origin' of the field cone 
ro = 3 x 10 9 cm is chosen to have the cone expansion of 
33.6°. The heliospheric expansion conserves the total 
number of electrons such that for scatter-free propaga- 
tion, J(r + r<j) 2 n(r)dr = const. 

The effect of the background electron density gradient 
on plasma waves is governed by the last term on the left 
hand side of Equation [51 Similarly to iKontarl (|2001al ) 
we define the characteristic scale of plasma inhomogene- 
ity, L = u} pe (dujp e /dr)~ 1 — 2n e (dn e /dr)~ 1 . This val- 
ues has to be larger than the wavelength of any plasma 
waves considered to remain within the Wcstzel-Kramcrs- 
Brillouin (WKB) approximation of geometrical optics. 

2.1. Initial conditions 

The electron distribution function is modelled using 
an instantaneous electron injection which is Gaussian 
in space with a characteristic size d. This electron 
distribution has a power-law spectrum in velocity, and 
hence in energy, as of ten observed in solar flares (e.g. 
iBrown fc Kontarl l2005!). f(v,r, t = 0) takes the form 



f(v,r,t = 0) =cxp - 



r 2 \ n b (2S-l) (u 

d 2 






Langmuir waves, and jl = V^^Jpe {v/vxe) exp ( — ^— j 



The electron beam is normalised to the electron num- 
ber density n\y. S represents the spectral index of the 
energy power-law and v m i n represents the minimum ve- 
locity used for the electron beam. 

The initial location of an electron beam (r — in the 
above equations) for the subsequent simulations is taken 
at a background plasma frequency of 500 MHz which 
corresponds to the height of 3 x 10 9 cm~ 3 above the 
photosphere. This is often interpreted as the typical fre- 
quency/locati on for an electron beam acceleration site 
in the corona (jAschwanden et aTlll995D . The spectral in- 
dex 5 was set to 3.5, corresponding to typical spectral in- 
dices above the brea k energy of in-situ m easured electron 
beams at the Earth ([Krucker et a l 2009) . The beam size 
was taken to be d = 10 9 cm. Electron thermal velocity 
was taken to be VTe = 5.5 x 10 8 cm/s, which corresponds 
to Maxwellian plasma with a temperature of 1 MK. The 
beam velocities will range between 3.6«Te ~2x 10 9 cm/s 
and 2 x 10 10 cm/s. Above the maximum velocity rela- 
tivistic effects become important. Plasma waves created 
near thermal velocity are absorbed by the background 
Maxwellian through Landau damping so i.&VTe is an ac- 
ceptable lower limit. 

The initial electron beam density is taken to be 1.1 x 
10 5 cm -3 which, together with S — 3.5, gives the to- 
tal number of electrons above 50 keV of 1.2(y / 7rd) 3 w 
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7 x 10 27 . This is a relatively small event in relation to ob- 
serve d number of electrons above 50 keV (jKrucker et al.l 
120071) . The instantaneous injection of the electron beam 
restricts the total injected electrons to small event sizes 
to keep the flux of electrons around 100 keV near the 
Earth in line with typica l values observed at 1 AU 
(|Krucker et al.|[2007l I2009T) . If we consider similar num- 
ber densities at the peak of a temporal Gaussian injection 
of order 10 3 s, the total number of electrons rises to 1 31 
in agreement with observations (jKrucker et alJl2007t ). 

The initial spectral energy density of the plasma waves 
is assumed to be at the thermal level 
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W(v,r,t = 0) = 



k B T e Upejr)' 

4vT 2 V 2 



lofi 



V Te 



(4) 



where T is the background plasma temperature, k B is 
Boltzmann constant and vx e is the background electron 
thermal velocity. This thermal level is formed setting 
dW/dt = for Maxwellian distribution of electrons with 
temperature T and ignoring electron collisions in Equa- 
tion (0). 

2.2. Heliospheric plasma density 

The background heliospheric plasma is modelled as 
a continuously decreasing background electron density. 
The background density model is found using the equa- 
tions for a s tationary sph erical symmetric solution for the 
solar wind (|Parkerlll958|) with a constant fo und by satel- 
lite m easu rements near the Ear th's orbit (|Mann et al.l 
[19991) . See lKontar fc Reidl (J2009I) for details. The model 
is static in time because the characteristic electron beam 
velocities are much larger than solar wind speeds. 

3. ELECTRON TRANSPORT THROUGH PLASMA 
WITH DECREASING DENSITY 

The initial electron distribution injected into the simu- 
lation is stable at t = but once the electrons are allowed 
to propagate through space, the distribution quickly be- 
comes unstable (df/dv > 0) to plasma wave genera- 
tion. This is related to the 'time of flight' effect first 
introduced bv iFilbert fc Kelloggl (| 19791 ). As the growth 
rate of plasma waves is velocity dependent, the initial 
power law distribution causes quasilinear relaxation to 
be important up to a certain velocity or corresponding 
break energy. Above this energy electrons are too dilute 
to generate any plasma waves and travel scatter free. 
Below this energy plasma waves are generated, relaxing 
the distribution function to a plateau in velocity space 
(df/dv rj 0) as energy from the electrons is transferred 
to the generated plasma waves. The in s tability forms 
a bea m-plasma structure (|Mel'Niklll995t iKontar' et al.1 
Il998| ) , between the electron beam and the corresponding 
induced plasma waves seen in Figure [TJ We note that the 
presence of non-thermal particles leads to spontaneous 
generation of waves even when df/dv < 0. However, 
the level of spontaneously generated waves is low and in- 
sufficient to change the electron distribution function at 
substantially. 

The electron beam travelling through an unperturbed, 
decreasing background plasma experiences a slowly vary- 
ing small density gradient. This decreasing gradient 
causes generated plasma waves to drift to smaller phase 
velocities as L is strictly negative. The wave energy at 
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Fig. 1. — Colour coded plot of the electron flux [cm 2 eV s] — 1 and 
spectral energy density (normalised by thermal level W (v, X, t = 
0)) of plasma waves for two moments of time. Distance and velocity 
are normalised by solar radii and thermal velocity respectively. The 
background plasma density is unperturbed. 



one point in space is thus redistributed over a wider range 
of phase velocities. Therefore some generated waves are 
taken out of resonance with the energetic electrons which 
induced them. Energetic electrons arriving later in time 
to this point in space are unable to reabsorb all the 
energy previously converted to plasma waves. Conse- 
quently over time, the total energy in the electron beam 
is decreased resulting in a deceleration of electrons below 
the maximum energy at which plasma waves are induced. 
As it was shown by IKontar fc Reidl (|2009f ). this effect 
leads to the formation o f a broken power-law in fluenc e 
often detected at 1 AU (JLinl 119851: IKrucker et all 1200 71) . 
The waves shifted towards lower velocities are eventually 
absorbed by the background thermal plasma via Lan- 
dau damping. The recent survey of in- situ measured 
impu lsive solar energetic electron events (jKrucker et all 
120091) suggests the break energies generally appear in the 
deca-keV range , in li ne with numerical simulations by 
IKontar fc Reidl (f2009h . 

4. ELECTRON TRANSPORT THROUGH PLASMA 
WITH DENSITY PERTURBED BY A SINE- WAVE 

The background electron density is a simplified model 
where only large scale radial expansion is taken into ac- 
count. The real inner hcliosphcre electron density has 
density fluctuations present at various smaller scales. To 
initially explore density fluctuations, a simple perturba- 
tion of the background plasma is added to the previous 
heliospheric density model. This perturbation takes the 
form of a sinusoid giving a new background density 



n(r) = no(r)[l + asin(27rr/A)] 



(5) 
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where a and A are the amplitude and wavelength of the 
perturbation respectively and rig is the original unper- 
turbed density. The initial value of the amplitude a is 
taken as 10 -2 while the wavelength A is taken as 10 10 cm. 
These values create a perturba tion which is within rea- 
sonable solar wind parameters (jCelnikier et alJ fl983). 

4.1. Distributions close to the Sun 

Close to the Sun, the radial drop of density is very 
sharp and plays the dominant role in density change. 
The small scale fluctuations are thus unable to gener- 
ate any positive density gradients. The drift of waves 
in velocity space is always to lower phase velocities 
which can be observed at the earlier time interval t = 
25 s (Figure [2]). The density fluctuations cause an 
increase or decrease in this drift of plasma waves to 
lower phase velocities. As the growth rate of plasma 
waves depends linearly upon the magnitude of plasma 
waves at any point in phase space, if the plasma in- 
homogeneity is too large then waves are shifted too 
fast and plasma wave production is suppressed (in line 
withlSmith fc Simdll97JlMuschietti et al.|[T985l : iKontarl 
I2001M ILedenev et alJl2004t iLi et al.ll2006bn . 

To compare the background plasma inhomogeneity 
with the level of plasma waves in any spatial location 
we consider the magnitude of wave energy density, found 
by 



Time 



25.0 seconds 



E w (r,t) = 



Wdk = 



W(r,v,t) 



dv. (6) 



The plasma wave energy density, E w (r, t) , close to the 
Sun at time t = 25 s is displayed in Figure [3] with the 
corresponding scale of the background plasma inhomo- 
geneity. The unperturbed case has been over plotted 
for comparison. Lines have been drawn to indicate the 
10 10 cm wavelength of sinusoid perturbation to the back- 
ground plasma. Periodic oscillation of the background 
plasma is evident together with the corresponding pe- 
riodic nature of the plasma wave energy density. The 
magnitude of E w (r, t) in the unperturbed case is gener- 
ally larger than the perturbed case, showing clearly the 
reduction in wave growth when the background plasma is 
significantly perturbed. As we get further away from the 
Sun (5R S compared with 2R S ) the radial drop of density 
plays a less dominant role allowing small scale fluctua- 
tions to become more important, seen in \L\~ . With 
this increased role, the small scale fluctuations increase 
the suppression of induced plasma wave energy density 
with respect to the unperturbed case. 

Despite fluctuations suppressing plasma waves, the 
perturbed case displays plasma wave energy density 
greater than the unperturbed case at peaks in its oscilla- 
tion. The instability of the electron beam which induces 
the plasma waves (df '/dv > 0) is not fully relaxed to 
thermal velocities in areas of space where plasma wave 
production is suppressed. Another striking feature of 
Figure [3] is the double peak and trough behaviour of 
E w (r, t) within one wavelength of background plasma 
fluctuation. 

The distribution of E w (r,t) in space is substantially 
different at the latter time of t = 100 s, shown in Figure 
[3j There is a larger discrepancy in magnitude between 
the unperturbed and perturbed case. Moreover, the sec- 
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Fig. 2. — Colour coded plot of the electron flux [cm 2 eV s] — 1 and 
spectral energy density (normalised by thermal level W(v, x,t = 
0)) of plasma waves for two moments of time. Distance and velocity 
are normalised by solar radii and thermal velocity respectively. The 
background plasma density has been perturbed with a sine wave. 



ond peak of E w (r, t) within one wavelength clearly seen 
at t = 25 s is suppressed at the later time of t = 100 s. 
The one remaining pronounced peak does not stay co- 
spatially with the small scale fluctuation wavelength but 
shifts backwards with respect to increasing distance from 
the Sun for this single point in time. Density fluctu- 
ations at distances rj 7R s become influential enough 
over the radial density decrease to generate some pos- 
itive background density gradients. A positive gradient 
causes plasma waves to move to higher phase velocities 
and causes the streaking seen at t = 100 s in Figure [5J 
Despite the plasma wave distribution being substantially 
different, the electron flux remains al most un c hanged , 
agreeing with the numerical results from lKontarl (|2001b[) . 

4.2. The role of Langmuir wave group velocity 

The group velocity of plasma waves, Zv\Jv, is small in 
magnitude, within the range 4x 10 7 cm/s to 4x 10 8 cm/s. 
At t = 25 s (Figure [3]) the removal of the group velocity 
term has minimal effect. Waves are moved in space by 
a small distance dependent upon the magnitude of the 
group velocity. The slower energetic electrons at the back 
of the beam produce waves with higer group velocity and 
hence the wave energy density is displaced further. 

At the later time of t = 100 s, E w (r, t) is substantially 
different when the group velocity term is not present, 
seen in Figure [3J There is a clear double peak and 
trough behaviour within one background density fluc- 
tuation wavelength. Without any group velocity, waves 
are unable to travel from points in space where the back- 
ground density structure favours wave growth to points 
where wave growth is suppressed. The simulation with 
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Fig. 3. — The plasma wave energy density E w (r) at two different 
times for background plasma which is unperturbed (purple), per- 
turbed (red) and perturbed without implementing group velocity 
(green). The corresponding magnitude of plasma inhomogeneity 
|L| _1 for unperturbed (black) and perturbed (blue) is plotted for 
comparison. The light blue diamonds are where the plasma inho- 
mogeneity is positive in magnitude. 

no group velocity also has a higher magnitude of wave 
energy density at its peaks than both the other simula- 
tions. 

The group velocity of plasma waves, despite being 
small, acts to move wave energy from points in space 
where plasma waves are strongly induced to points in 
space where they are suppressed. This has a spatial 
smoothing effect on the induced plasma wave energy den- 
sity. 

4.3. Amplitude of fluctuations 

The amplitude a of the density fluctuations directly 
varies the background electron plasma density. The mag- 
nitude of this factor near the Earth can be found from 
observational results. It has been measu red using the 
Inter national Sun-Earth Explorer (ISEE) (jHarvev et al.1 
11978( 1 propagation experime nt with scintillation tech- 
niques (jCelnikier et al .1119831 ) that the background elec- 
tron plasma density near the Earth varies by about 10%. 
This would give the amplitude of a < 0.1. Therefore we 
consider a in the range 10~ 3 < a < IO -1 . The wave- 
length of the perturbation was taken as A = IO 10 cm. 

As a increases in magnitude, the oscillation in wave 
energy density increases. Similarly as a decreases in 
magnitude, the oscillations in wave energy density de- 
crease such that as a — > the wave energy density tends 
to the state where no perturbations arc present in the 
background electron plasma density. This can be seen 
in Figure E] in the plasma inhomogeneity, |i| _1 . As a 
decreases to IO -3 , the plasma inhomogeneity does not 
vary as much and L remains negative. 



The variation of a in Figure U shows how the mag- 
nitude of the plasma inhomogeneity affects wave gener- 
ation. If the fluctuations are too large, plasma waves 
drift in phase velocity too fast from the beam and are 
unable to build up. This suppression can clearly be seen 
when a = 10 . Most spatial areas have large values 
of |L| _1 and corresponding low values of wave energy 
density. Conversely, when a = 10~ 3 , the small scale 
fluctuations are small and wave energy density is able to 
build up to high magnitu des. This suppression agrees 
with previous theoretical (|Mehosd |1982t iMelrose et al.l 
119861) and numerical wor k on Langmuir wave generation 
in non-uninform plasma (jKontarll2001b[) . 

4.4. Wavelength of perturbations 

The wavelength of density fluctuations A has a strong 
effect on the local scale of plasma inhomogeneity, L, 
through dn/dr having one term inversely proportional 
to A. Density fluctuations have been measured at a va- 
riety of different length scales from IO 12 cm down to 
IO 6 cm (iNeugebauer et~all Il978t iCelnikier etHI 119871 : 
iKellogg et al.l I2009T T We have varied A in the range 
IO 9 cm < A < IO 11 cm which is close to the r ange of 
fluctuations presented by ICelnikier et al.l (|1987| ). The 
amplitude was set to a = 10~ 2 , similar to the previous 
section for comparison reasons. 

As A increases in magnitude, the oscillation in den- 
sity inhomogeneity decreases such that as A — > oo, the 
wave energy density tends to the state where no pertur- 
bations are present in the background electron plasma 
density. This can be seen from Figure [4] in the case 
where A = 10 n cm and the density inhomogeneity is 
very smooth. Conversely, as A decreases, the magnitude 
of L^ 1 increases while the sign of L^ 1 fluctuates rapidly. 

We can see from Figure 0] that when A is large, the 
induced plasma wave energy density resembles the un- 
perturbed case. When A is small, the large magnitude 
of L^ 1 causes waves to shift in velocity space faster. At 
any spatial point waves are present with a far greater 
range of phase velocities, however, their magnitude is 
much decreased. This means there exists a decreased 
level of plasma waves at points in phase space where the 
electron beam is present. The growth factor of plasma 
waves, responsible in the kinetic equations for converting 
electron beam energy to plasma wave energy, is propor- 
tional to the level of plasma waves. The decreased level 
of plasma waves in areas of phase space where the elec- 
tron beam is present causes less energy to be transferred 
from electron beam to plasma wave and is the reason 
for the smaller wave energy density observed in Figure [4] 
when A = 10 9 cm. 

5. ELECTRON TRANSPORT THROUGH PLASMA 
WITH POWER-LAW DENSITY FLUCTUATIONS 

The power spectrum of density fluctuations observed 
in the solar wind density follows a simple, Kolmogorov 
type power law near the Earth with spectral index 
near to 5/3. A similar spectrum index of perturba- 
tions h as been observed both w ith scintillation tech- 
niques (ICelnikier et al.|[T983lll987T) and with in-situ mea- 
surem ents (jNeugebauer et al.l 119781 : IKellogg fc Horburvl 
120051 ). The spectrum has been observed to steepen at 
small wavenumbers around 10 8 cm. To model small scale 
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Fig. 4. — Top: Energy density of Langmuir waves for a = 10 1 
(green), IO" 2 (red), 10~ 3 (blue). A = 10 10 cm. Bottom: Wave 
energy denisty for A = 10 11 cm(green), 10 10 cm(red), 10 s cm(blue). 
a = 10~ 2 . The background plasma inhomogcneity L(r) for each 
simulation in the appropriate colour is shown in lower panels. 



density fluctuations many perturbations of the back- 
ground plasma are introduced, so the density is 



n(r) = n (r) 



N 



1 + C J2 Afsm(2^/A n + 4> n ) 



(7) 



for N perturbations where no(r) is the original unper- 
turbed density. A n arc the wavelengths of density pertur- 
bations with <p n as their random phase. C is a constant 
which normalises the density fluctuations given by 



C = 



2(An(r 



N 



WE^xA: 



(8) 



where (n(r)) is the mean density. The root mean squared 
deviation of the density , y/(An(r) 2 ), from observational 



values near the Earth was taken t o be 0.4 c m" 3 or 10% of 

the mean density. The quantity J /JJlVyJ , the fractional 

density fluctuations, is a measure of the turbulent inten- 
sity of the background plasma. From Equation [8] this 
quantity is radially independent giving a constant tur- 
bulent intensity from the Sun to the Earth. We can then 
model the radial variation of turbulent intensity with 



' (An(r)2) f n (lAU) \ 
(n(r))* { n (r) J 



(An(r 1AU )*) 

{n{r 1AU )) 2 



(9) 



where ip > determines the degree at which the den- 
sity fluctuations become less dominant with ip = cor- 
respoding to no radial variation. For simplicity, we will 
reference the fractional density fluctuations as An/n. We 
considered the range on A to be IO 7 cm < A < IO 10 cm 
which is within the inertial range of solar wind turbu- 
lence. Larger values of A have a minor effect and the 
amplitude of waves shorter than A ss IO 7 cm is small. 
The random phases < </> < 2n ensure the amplitudes 
of density fluctuations have a Gaussian distribution. 

A constant level of An/n thoughout the inner helio- 
sphcre is found by setting ijj = 0. We set An/n = 10% 
which is within the observ ed range of values near the 
Earth (|Celnikier et al.|[l987f ). Figure [5] shows the density 
inhomogcneity and corresponding plasma wave energy 
density close to the Sun. The high level of inhomogcne- 
ity caused by the small scale fluctuations greatly sup- 
presses plasma wave spatial build-up compared to the 
unperturbed case. This suppression of plasma waves for 
An/n — 10% close to the Sun can prevent the high level 
of plasma waves required for type III solar radio emis- 
sion. 

To vary the level of fluctuations from the Sun to the 
Earth, we set ip > 0. Figure [5] shows the density inho- 
mogeneity and corresponding plasma wave energy den- 
sity close to the Sun with An/n = 10% at 1 AU and 
ijj = 0.5. Comparing the plasma wave energy den- 
sity with the earlier simulations which assumed constant 
An/n, we can see a much larger magnitude of plasma 
wave energy density being induced from the electron 
beam. The reduced contribution from the small scale 
fluctuations allows build up of plasma wave energy den- 
sity. This high level of wave energy density is required to 
see the recorded brightness temperatures associated with 
type III radio emission. Moreover, the spatial spread of 
plasma waves is much less sporadic than the produced 
level of wave energy density in the simulation with high 
level of fluctuations. The electron beam and plasma 
wave distribution can be seen in Figure [6] Plasma waves 
are no longer perturbed in a periodic fashion but are 
pseudo-random in space. The pseudo-random nature of 
the spikes in plasma wave energy density leads to similar 
clumpy behavio ur of plasma waves observ ed in-situ by 
spacecrafts (e.g. lGurnett fc Andersor]|1976l ). 

6. ELECTRON SPECTRA NEAR THE EARTH 

Previous work (jKontar fc Reidl I2009T) has shown the 
generation and absorption of plasma waves coupled with 
the effect of the background plasma inhomogcneity can 
change the electron beam energy distribution. A broken 
power-law in fluence spectrum can be formed from an ini- 
tially single power-law distribution. The break at which 
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Fig. 5. — The plasma wave energy and corresponding plasma in- 
homogeneity when density fluctuations have a power law spectra 
in frequency space and An/n = 10% at the Earth Top: The fluc- 
tuations arc constant from the Sun to the Earth (tp = 0). Bottom: 
The fluctuations increase from the Sun to the Earth (ip = 0.5). 
Both graphs arc over plotted with the unperturbed case (green). 
The plasma inhomogncity is plotted for unperturbed case (black) 
and perturbed case (blue) with light blue diamonds for positive 
values. 
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Fig. 6. — Colour coded plot of the electron flux [cm 2 eV s] -1 and 
spectral energy density (normalised by thermal level W(v,x,t = 
0)) of plasma waves. Distance and velocity are normalised by solar 
radii and thermal velocity respectively. The small scale fluctuations 
have a power law spectra in frequency space where the fluctuations 
increase from the Sun to the Earth with ip = 0.5. 



the two power-laws connect is the maximum velocity the 
electrons were able to induce plasma waves through a 
resonant interaction with the background plasma. The 
spectrum below the break is flattened during transport 
because the electron beam is unable to re-absorb all the 
energy transferred to plasma waves due to background 
plasma density gradients. 

Introducing small scale density fluctuations into the 
background plasma changes its properties and should 
consequently change the spectrum of the electron distri- 
bution function. Whilst changes in the electron spectrum 
are not visible on short scales (a few relaxation times), 
the fluence spectrum at the Earth shows a noticeable de- 
pendence upon the level of fluctuations introduced to the 
simulation. Figure [7] shows the fluence spectrum of the 



electron beam at the Earth for five different amplitudes 
of fluctuation within the range 10 -3 < a < 10 — 1 . As 
shown earlier, the small scale density fluctuations sup- 
press the generation of plasma waves. This decreases 
the amount of energy transferred through resonant in- 
teraction from the electron beam to the plasma waves. 
With less total energy, a smaller amount of energy in 
plasma wave form can drift to higher or lower phase ve- 
locities and not be re-absorbed by the electron beam. 
The amount of deceleration the electron beam can un- 
dergo due to plasma waves drifting is decreased, causing 
a reduction in the flattening of the fluence electron spec- 
trum. This means when a is larger, the fluence spectrum 
below the break energy has a higher spectral index (Fig- 
ure [8]). Similar behaviour is demonstrated by the fluence 
spectrum of the electron beam at the Earth for four dif- 
ferent wavelengths of small scale fluctuations within the 
range 10 8 cm < A < 10 11 cm, shown in Figure [7J The 
same lack of wave generation for small A reduces the 
deceleration the electron beam undergoes and hence re- 
duces the flattening of the fluence spectrum (Figure [8]). 

Despite the change in spectrum near the Earth, the 
electron distribution function does not share the same 
sensitivity to the structure of the background electron 
density as the plasma waves (See Figures [T] and [2]) . The 
simulation with perturbed plasma does however show 
small changes, most noticeably in the tail of the electron 
distribution. A positive spatial gradient in background 
plasma causes plasma waves to drift to higher phase ve- 
locities. This drifting of waves in velocity space allows 
their energy to be re-absorbed by the tail of the beam 
such that electrons are accelerated to higher energies. It 
is the opposite effect of the negative density gradient tak- 
ing plasma wave energy away from electrons and forming 
a broken power-law near the Earth. This acceleration of 
electrons causes the noticeable bump around 10 — 20 kcV 
in Figure seen for simulations with higher spectral in- 
dices below the break energy. The bump becomes more 
prominent for small A, high a or more generally when 
the background density fluctuations are more effective 
at moving wave energy to higher phase velocities. 

Extending the density fluctuations to multi- wavelength 
model, a Kolmogorov type power-law is assumed where 
(An/n) 2 ~ A 5 ' 3 with An/n remaining radially constant 
(V' = 0). Figure [7] shows the fluence spectrum at the 
Earth for four different turbulent intensities 0.01% < 
An/n < 10%. The larger An/n is, the greater the sup- 
pression of plasma waves and hence the higher the spec- 
tral index below the break energy of the fluence spectra 
(Figure [5]). The signature bump can be seen in the spec- 
tra around 10 — 20 keV, again caused by the acceleration 
of electrons at the back of the beam due to plasma waves 
drifting to higher phase velocities. 

The electron beam fluence spectra for density fluctu- 
ations An/n changing with distance (reaching 10% at 
1 AU) are displayed in Figure [7J for four different val- 
ues of tp within the range < tp < 0.8. The decreased 
presence of fluctuations near the Sun (ip > 0) increases 
the amount of induced plasma waves which decreases the 
spectral index below the break energy, shown in Figure 
[5J For all simulated values of tp > no bump in the flu- 
ence spectra is present. A reduced level of fluctuations 
near the Sun decreases positive density gradients which 
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Fig. 7. — Fluence of the electron distribution function near 
the Earth. Top Left: Five simulations with a = 10 — 1 (black), 
lCT 1 ' 5 (purple), 1CT 2 (blue), 10" 2 ' 5 (green), and 1CT 3 (red). 
A = 10 10 cm. Top Right: Five simulations with A = 10 s cm 
(black), 10 cm (purple), 10 cm (blue), 10 cm (green) and un- 
perturbed (red), a = 10 -2 . Bottom Left: Four simulations for 
multi-scale fluctuations with An/n of 10 % (black), 1 % (blue), 
0.1 % (green) and 0.01 % (red) of the mean background density. 
Bottom Right: Four simulations for multi-scale fluctuations which 
decrease in power close to the Sun for ip of (black), 0.3 (blue), 
0.5 (green), 0.8 (red). 




FlG. 8. — The spectral index of a power law fit between 4 and 40 
kcV for the fluence spectra of electrons near the Earth. Top Left: 
Spectral index versus the amplitude of density fluctuation. Top 
Right: Spectral index versus the wavelengths of density fluctuation. 
Bottom Left: Spectral index versus multi-scale level of fluctuations. 
Bottom Right: Spectral index versus ip, the radial degree at which 
density fluctuations become less dominant. 

subsequently decreases the acceleration of electrons from 
re-absorption of plasma waves. A smoother increase in 
spectral index below the break energy for increasing ip 
can thus be seen in Figure [5] 

7. DISCUSSION AND CONCLUSIONS 

The simulations show that fine structure of the back- 
ground solar wind electron density caused plasma waves 
to be suppressed, with larger amplitudes and smaller 
length scales of density fluctuations having the largest ef- 
fect. This increased suppression for larger amplitudes is 
similarly observed for higher levels of turbulence (An/n) 



with Kolmogorov type density fluctuations. 

For high levels of turbulence near the Sun, An/n = 
10%, wave production by the electron beam became 
no longer sufficient for the generation of type III radio 
bursts. It is possible to induce more plasma waves by 
increasing the density of the electron beam. This solu- 
tion requires at least two orders of magnitude more elec- 
trons, causing the beam to have around 1% of the density 
of the background plasma. Such high density electron 
beams become problematic when considering simultane- 
ous HXR bursts assuming the upward electron beam has 
only 0.2% of the downward electron b eam density, found 
above 50 kcV in lKrucker et al.l (|2007| ). 

Increasing the level of plasma waves near the Sun with- 
out increasing beam density, the amplitude of density 
fluctuations can be reduced. We implemented a radial 
dependence with closer distances to the Sun have a de- 
creased turbulent intens ity. This is seen in observat ional 
scintillation techniques tiWoo et al.lll995HWoolll996t) and 
Helios data (|Marsch fc Tulll990l) in the fast solar wind. 
The observed val ues for An/n are as low as 0.3% at dis- 
tances < 0.1 AU (jWoo et al.lll995[ ). A much higher mag- 
nitude of plasma wave energy density was achieved close 
to the Sun with smaller levels of fluctuations. 

To estimate how density flucutations might radially 
evolve, we varied the initial conditions of the simula- 
tions. We used a variety of different initial electron 
beam spectral indices (8 in Equation [3]) and different 
radial dependence of density fluctuations (ip in Equation 
[5]). Using the resulting fluence spectra near the Earth 
for each simulation, we compared the spectral indices 
above and below the break energy (Figure IH]). The spec- 
tral index becomes smaller below the break energy for 
larger values of ip. We have also overplotted the cor- 
relation of spectral indices above and below the break 
energy of pe ak flux measurement s taken from a statis- 
tical survey dKrucker et al.l [20091 ) of impulsive electron 
events detected by the three-dimensional Plasma and 
Energetic Partic les experiment on the WIND spacecraft 
(jLin et al.lll995| ). A level of fluctuations with ip around 
0.25 would give a similar correlation to the observational 
line. We no te, however, that the observational line pre- 
sented from lKrucker et all (J2009D fitted a large scatter of 
data points. The ratio of low:high spectral index for all 
simulated results presented in figure [9] lies between 0.42 
and 0.58 which is w ithin the narrow range presented in 
iKruckeret all ([200l . 

A variety of simulation variables can affect the energy 
of the spectral break at the Earth: the model of radial 
background density decrease, the density fluctuations, 
the initial spectral index of the beam, the number density 
of injected electrons, the spatial distribution of injected 
electrons, the temporal nature of the injection, and the 
initial coronal background density where the electrons 
are injected. The spectral index below the break energy 
of the resultant double power-law in fluence spectra near 
the Earth is increased when density fluctuations have a 
larger effect on the level of induced plasma waves. It is 
important to note, however, the spectra below the break 
energy is only approximately a power-law. The presence 
of density fluctuations causes fine structure to be present. 
A bump around 10-20 keV was found, caused by accel- 
eration of electrons at the back of the beam through ab- 
sorbed plasma waves. The onset of this bump appears 
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Fig. 9. — Comparison between the high and low spectral index 
of fluence spectra of electrons near the Earth. The dashed purple 
line is t he best fit to the ob servational data of peak flux spectral 
indices IIKruckcr ct al.ll2u09t) . 

to be close to the Sun where plasma wave energy density 
is high. The magnitude of this bump is reflected in the 



size of the spectral index error bars in Figure |H] with a 
larger bump corresponding to a larger error. With the 
prospect of Solar Orbiter and Solar Probe Plus, it is very 
attractive to extend these studies further to understand 
the spectral evolution of the electron beam between the 
Sun and the Earth. 
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